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Detailed Protocol for scVelo for Single-nuclei RNA seq of the mouse placental labyrinth development 


1. Single cell velocity trajectory analysis is performed in python using the Velocyto package (LaManno, 2018). Install velocity according to the instructions 
from the developers (http://velocyto.org/velocyto.py/install/index.html). 


2. Generate intron/exon counts from expression matrices generated by 10x Genomics CellRanger3.0 using the run10x command in velocyto. This will 
generate loom files for each separate sample. 

3. Concatenate these matrices into one file for all timepoints using the loompy package loompy.combine function in python. This can be done in jupyter 
notebooks. 

4. To analyze the velocity for cells identified and clustered according to prior analyses performed in Seurat, convert the Seurat data to the loom format 
using the as.loom command in Seurat in R (https://satijalab.org/seurat/articles/conversion_vignette.html). 

5. Install scvelo according to the instructions from the developers (https://scvelo.readthedocs .io/installation.htm!). 

6. The scVelo analysis script is attached, and the steps are outlined below. 


7. To filter the output loom from velocyto for only the cells used in the trophoblast analysis from Seurat, use the utils. merge function within scVelo. For our 
data, the concatenated output from velocity (AllStages_velocyto_output.loom) and the integrated trophoblast dataset 
(AllStages_TrophoblastNuclei.loom) are available for download here: https://figshare.com/projects/Single_nuclei_RNA- 
seq _ of mouse placental labyrinth _development/92354. 


8. Use the scvelo.pp.moments function to calculate the moments for velocity estimation, using the default settings of using 30 PCs and 30 neighbors. 
Then the proportion of spliced, unspliced, ambiguous reads was determined using scvelo.utils.show_proportions. 
9. The velocity was then calculated using scvelo.tl.velocity command, specifying the mode=’stochastic’. The velocity graph was then computed using 
scv.tl.velocity_graph with default parameters. 
10. The data were then split by timepoint (E9.5, E10.5, E12.5, E14.5) and the velocity embeddings visualized using the scvelo.pl.velocity_embedding 
function. 
11. The magnitude of the velocity vectors was estimated from the velocity embeddings for each cell given by the ‘velocity_umap_cell_embeddings’, by | 
calculating the distance from the cell in umap space from the starting coordinate or each cell (0,0) to the coordinate of the velocity embedding provided Ee 
by the ‘velocity_umap_cell_embeddings’. The units for this measurement are arbitrary and specific to the UMAP space. 
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